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Abstract 



We calculate the cb mass spectrum, the splitting values and some other 
properties in the framework of the semi-relativistic equation by applying the 
shifted large-N expansion technique. We use seven different central potentials 
together with an improved QCD-motivated interquark potentials calculated to 
two loops in the modified minimal-subtraction (MS') scheme. The parameters 
of these potentials are fitted to generate the semi-relativistic bound states of 
c6 quarkonium system in close conformity with the experimental and the 
present available calculated center-of-gravity (c.o.g.) data. Calculations of 
the energy bound states are carried out up to third order. Our results are in 
excellent fit with the results of the other works. 



I. INTRODUCTION 

The spectroscopy of the c6 system have already been widely studied before in the framework 
of the heavy quarkonium theory [1]. After discovery of the was reported in 1998 by 
the (CDF) collaboration [2], the observed mass Mb, = GAO^^il GeV has inspired new 
theoretical interest [1,3-9]. Bound state masses have been estimated for the system 
which consists of heavy quarks [1,3-9]. Hence, it can be reliably described by the use of the 
methods developed for the cc and the bb spectra. 

Quite recently, the revised analysis of the Be spectroscopy has been performed in the 
framework of the potential approach [1,4,7,9] and QCD sum rule [3,6]. Kwong and Rosner 
[7] predicted the masses of the lowest cb vector (triplet) and pseudoscalar (singlet) states 
using an empirical mass formula and a logarithmic potential. Eichten and Quigg [1] gave 
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a more comprehensive account of the energies and properties of the c6 system that was 
based on the QCD- motivated potential of Buchmiiller and Tye [10]. Gershtein et al. [8] 
also published a detailed account of the energies and decays of the cb system using a QCD 
sum-rule calculation. Baldicchi and Prosperi have computed the c& spectrum based on an 
effective mass operator with full relativistic kinematics [6]. Moreover, they have fitted the 
entire quarkonium spectrum. Fulcher et al. [4] extended the treatment of the spin-dependent 
potentials to the full radiative one- loop level and thus included effects of the running coupling 
constant in these potentials. He also used the renormalization scheme developed by Gupta 
and Radford [11]. Very recently, we have applied the shifted large-N expansion technique 
(SLNET) [12-14] to study the c6 system in the context of Schrodinger equation [9]. 

In this work, we choose a group of seven central potentials [1,13-24] with a strong coupling 
constant to fit the spectroscopy of the present spin-averaged energies of the existing 
quarkonium systems. We insist upon strict flavor-independence of their parameters. We 
also extend this study to an improved QCD-motivated interquark potential calculated to two 
loops in the modified minimal-subtraction (M^*) scheme [25,26]. Such a potential contains 
an expansion to ag to at least 2-loop order and use of nontrivial methods in the interpolation 
between long-distance (string constant) and short-distance behaviors [10]. 

Since one would expect the average values of the momentum transfer in the various 
quark- antiquark states to be different, some variation in the values of the strong coupling 
constant and the normalization scale in the spin-dependent should be expected. 

In order to minimize the role of flavor-dependence, we use the same values for the coupling 
constant and the normalization scale for each of the levels in a given system and require 
that these values be consistent with a universal QCD scale. 

The contents of this article will be as follow: In section II, we present the solution of the 
spinless Salpeter equation using the SLNET for the c6 Salpeter spin-averaged binding en- 
ergies. Section III presents the proposed phenomenological central and the QCD-motivated 
potentials. Summary and conclusions are presented in section IV. 
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II. WAVE EQUATION 



The relativistic wave Salpeter equation [14] is constructed by considering the kinetic energies 
of the constituents and the interaction potential. The spinless Salpeter (SS) for the case 
of two particles with unequal masses rUg and tuq, interacting via a spherically symmetric 
potential V{r) in the center-of-momentum system of the two particles is given by 

sJ-^N + ml + V{r)-M{qQ) 

where the kinetic terms involving the operation >/— Ajv + m| are nonlocal operators and 
'0(r) = Y^^jn{9, (f))Rn,i{r) denotes the Salpeter's wave function. For heavy quarks, the kinetic 
energy operators in Eq. (1) can be approximated, (cf. e.g., Jaczko and Durand of Ref. [27]), 
as 

^ + (2) 

i=q,Q ^ ' 

mm f m m \^/^ 

where a — ^''.'"'^ denotes the reduced mass and n — [ — '^'^'^^ \ . u is a useful mass 

^ mq+mq ' ymqmQ—Zn-^ J ^ 

parameter. This SS-type equation retains its relativistic kinematics and is suitable for 
describing the spin-averaged spectrum of two bound quarks of masses rUq and mq and total 
binding mass M{qQ). The Hamiltonian containing the relativistic corrections up to order 
(v^/c^) is called as the generalized Breit- Fermi Hamiltonian (cf. e.g., Lucha et al. of Ref. 
[21]). Therefore, the spinless Salpeter equation can be written (in units h = c = 1) [14] 

{"^ " + ^^"^} ^''^"^ ^ E^,iRnAr): (3) 
where E^^i — Mn,i{qQ) —mq — niQ refers to the Salpeter quark binding energy with Mn,i{qQ) 
is the semirelativistic-bound-state masses of the non-self-conjugate atomlike meson such as 
cb meson and A^r = V^.^ Furthermore, in order to obtain a Schrodinger-like equation, the 
perturbed term in Eq. (3) is treated using the reduced Schrodinger equation [29] 

^This approximation is correct to O^v^/c^). The term in (3) should be properly treated as a 
perturbation by using trial wavefunctions [28] . 
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(r) = 0, 



(1) 



(4) 



with — A%, and consequently one would reduce Eq. (3) to the Schrodinger-type form 

[14] 



{-^ - 1^ [Ell + V\r) - 2E^,V{r)] + V(r)| i?„,(r) = E^Mr). 

Further, the A^— dimensional space operator in the spherical polar coordinates is 

92 N-1 d L%_^ 



dr'' 



r dr 



with Llf_^ = I {I + N — 2) . After employing the following transformation 



Un,i[r) 

AN-l)/2 ■ 



and proposing the shifting parameter k — k — a with k — N + 21, one obtains [14] 



92 {k -l){k- 3) 



2 



and 



{k -l){k- 3) 92 {k -l){k-3) d 



Q^4 



+ 



dr 



(5) 



(6) 



(7) 



(8) 



{k - if {k - 3f - 24 (A; - 1) {k - 3) 
16r^ ■ 



(9) 



Thus, using Eqs (7) and (8) and after a lengthy manipulation but straightforward, we may 
write Eq.(5) in more simple and explicit form as 



1 , [k-{l-a)]\k-{3-a)] ^ ^^^^^^^^ Wr,,{rr 



2iidr^ 



2m 



Un,i{r) = 0, (10) 



with 



Wr^,i{r) = V{r)-Er,,i, 



(11) 



and the effective mass 
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(12) 



It is worthwhile to note that the expression (10) is in complete agreement with the expansion 
formula made by Durand and Durand in Ref. [27]. The perturbation term, W„,;(r)^; that 



condition is verified by the confining potentials used to describe heavy— quark systems except 
near the color— Coulomb singularity at the origin, and for r ^ oo. However, it is always 
satisfied on the average as stated by Durand and Durand [27]. We now proceed to solve 
Eq. (10) by applying the SLNET, obtaining an exact result in series for the energy [12-14]. 
This is a fundamental feature that allows one to attack problems that do not involve small 
coupling constant or Hamiltonians without a solvable strong term. Such an approximation 
has been used for the solution of nonrelativistic as well as relativistic wave equations with 
spherically symmetric potentials [12-14] yielding sufficiently accurate results. Imbo et al. 
[12] showed the high accuracy and the applicability of this method to a large number of 
spherically symmetric potentials. In the SLNET [9] it is convenient to shift the origin of 
coordinates to r = ro (or x — 0) by defining 



is, (w^/c^) term in Eq. (10) is significant only when it is small (i.e., Wn,i{r)/m' <^ 1). This 



X — 



(13) 



Expansions about this point in powers of x and k yield [9,13,14] 




(14) 




(15) 



and also 



1 



oo 




(16) 



Further, by substituting Eqs. 

\ I dp \ k 



(14) through (16) into Eq. (10), one gets 



4:11 dx"^ L16a* 
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Q dri ) 3\ m'Q 



h\ K ) 



with the eigenvalues 



Q 



k 



^0 



m 



I „ '^EnEo E? \ 1 
+ ( ^2 + + ^ T 

m m J k 



(17) 



+ ( + + ) p + 



m 



(18) 



The parameter Q is an arbitrary scale, but is to be set equal to k to rescale potential 
after the problem is solved formally in the 1/k expansion. Thus comparing Eqs. (17) and 
(18) with its counterpart Schrodinger— like equation for the one— dimensional anharmonic 
oscillator problem which has been investigated in detail for spherically symmetric case by 
Imbo et al. [12] . The final analytic expression in the 1/k expansion of the energy eigenvalues 
appropriate to the SS particle is 



S 



k 



16// 



Q 



m'Q 



m'Q 



+ 



1 



m'Q 



1 



2rlEsV{ro) 
m'Q 



+ 



1 



(19) 



but Ur is to be set equal to 0, 1, 2, • • • . The quantities /S^^^ and (3^'^^ appearing in the correction 
to the leading order of the energy expression are defined and listed in the Appendix A. 

Thus, comparing the terms of Eq. (18) with their counterparts in Eq. (19) and equating 
terms of same order in k yields the leading energy [14] 



where m^u = ///m'and ro is chosen to minimize the leading energy, Eq. (20), that is, [9,13-14] 

(21) 



dro 



drl 



Therefore, Tq satisfies the following expression 



rtV\T,)[— + 



m 



Q 



4 ' 16r^m^ 



1/2 



Q 



16m„ 



(22) 



Further, to solve for the shifting parameter a, the next contribution to the energy eigenvalue 
is chosen to vanish, (i.e., Ei — 0),(cf. e.g. Refs. [9,12-14]), which gives 



a = 2 - 4/x(l + 2nr)u;, 



(23) 



where cu is given by 



= ^ [3 + roV"{ro)/V'{ro) - IQr'.m.V'irof/Q] 
and Q in Eq. (22) can be rewritten in a more convenient form as 

Q = 8m4roV(ro)]'(l + A), 



(24) 



(25) 



with 



A = Wl + 



roV^'(ro) 



(26) 



Therefore, solving Eqs. (23) through (26) together with Q — k ^ yields 
1 + 2Z + A^l{2nr + l)uj = 2rlV\ro) (2m^ + 2m^Xf'^ , 



(27) 



which is an explicit equation in tq. Once ro is determined via Eq. (27), Eq can be obtained 
via Eq. (20), E2 and £'3 are also obtained by solving Eqs. (18) and (19). Consequently, the 

general expression for the quark binding energy is 



End — -En + 



1 



(28) 



which is an elegant algebraic expression that gives a rapidly convergent binding energy value 
with high accuracy. Now, in the framework of our semi-relativistic independent particle 



model the mass levels for an atomlike qQ meson where a light quark q is moving around an 
almost fixed heavy anti-quark Q, the Salpeter bound-state mass can also be easily obtained 
from the expression of quark binding energy (28) as 

MnMQ) =mq + mQ + 2En,i. (29) 

Here, of course, we do not consider the recoil effects of the heavy anti-quark Q. Finally, 
making use of Eqs. (27) through (29), one can resolve Eq. (10) for different types of central 
potentials taking into account the spin dependent terms Vsd as correction terms to the static 
potential. We also consider a common parameter set for the upsilon and charmonium mass 
spectra in a flavour-independent case. 



III. SOME POTENTIAL MODELS 



A. Static potentials 



The cb system that we investigate in the context of SS wave equation is often considered 
as nonrelativistic system, and consequently our treatment is based upon Salpeter equation 
with a Hamiltonian 

H^Ho + H, + Vsd, (30) 

where Hq — — Ajv/2/x-|- \^(r) denotes the Schrodinger Hamiltonian, Hi = —A%/8r)^ denotes 
the perturbation term and Vsd is the spin-dependent term [1,4,9,21] given by^ 



^To the moment, the only measured splitting of nS'-levels is that of rjc and J/i^, which allows 

us to evaluate the so-called (c.o.g.) data using M^(IS) = (3Mj/^ + Mrjc)/^ and also M(nS) = 
My(nS) - {Mj/^ - Mr,c)/4n [15,21]. 



Spin-independent relativistic corrections are not included. For the purpose of making some 
preliminary estimates of the energies of the lowest two S-states of the B^. system, it is nec- 
essary to consider only the spin-spin part of the spin- dependent potential since these are 
S-states. Our solutions to the SS equation are generated numerically for the central poten- 
tials. The effects of the spin-dependent parts are added as a perturbation as improvement 
for splitting the lowest S-state of the spin-averaged energies. The potential parameters in 
this section are all strictly flavor- independent. The potential parameters and the constituent 
mass parameters are fitted to the low-lying energy levels of charmonium and upsilon systems. 
So, like most authors (cf. e.g. [1,4,8,9,21]), we determine the coupling constant as{mc) from 
the well measured experimental charmonium hyperfine splitting of the IS(cc) state value 
of AE'hf — Mj/^ — M^^ = 117 ± 2 MeV. The numerical value of ag is dependent on the 
potential form and found to be compatible to the other measurements [1,3-4,6-9,15]. The 
observed Mj/^ — M^^ hyperfine splitting fixes as for each potential. The perturbative part 
of such a quantity was evaluated at the lowest order in a^. Baldicchi and Prosperi [6] used 
the standard running QCD coupling expression 



with Uf — 4 and A — 200 MeV cut at a maximum value q;s(0) — 0.35 to get the right 
Mj/^ — Mj^^ sphtting and to treat properly the infrared region [6] . Whereas Brambilla and 
Vairo [3] took in their perturbative analysis 0.26 < ^^(/x = 2 GeV) < 0.30. 
We consider now a class of static potentials of general form 



had been previously proposed by Lichtenberg et al. [17] in which the additive constant Cq 
may be of either sign. It was necessary to add a substantial flavor-independent constant to 
each potential for each quark-antiquark system. The constant Co in Eq. (33) is determined 
from the fit to the spin-averaged data (SAD) or center-of-gravity (c.o.g.) mass of the IS 
state (i.e., M(1S) = 3068 MeV) from the PDG [30]. Since the model is spin independent 





(32) 



V{r) = -^r"^ + + Co; < 7 < 1, A, 5 > 



(33) 
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and as the energies of the singlet states of quarkonium famihes have not been measured, a 
theoretical estimate of these unknown levels introduces uncertainty into the calculated SAD 
[9,31]. The best parameter fittings to (c.o.g.) data are obtained via the chi-square value or 
the least-square fit x ' ^j^" — , where AMj refer to the errors in the experimental 

i ' - 

data. The static quarkonium potentials are monotone nondecreasing, and concave functions 
which satisfy the condition 

V\r) > and V"{r) < 0. (34) 
This class of potentials, in the generality (33), comprises the following family of potentials: 



1. Cornell potential 

The QCD-motivated Coulomb-plus-linear potential (Cornell potential) has the form [24] 

rv 

Vc{r)^-^ + Kr + Co, (35) 
r 

The main drawback of this potential is that the cc and bb states lie in an intermediate 
region of quark separation where neither limiting forms of (34) should be valid. 

a. Coulomb and linear potential of Eichten et al. [1,22,24]: with etc — 0.52, n — 
0.1756 GeV^ Co = -0.8138 GeV, = 1.840 GeV,and 

mfe = 5.232 GeV. We label this potential by Cornell 1. 

b. Coulomb and linear potential (35) of Hagiwara et al. [18,22]: with etc = 0.47, k ~ 
0.19 GeV^ Co = 51 MeV, rric = 1.32 GeV and = 4.729 GeV. We label this potential by 
Cornell 2. 

2. Song-Lin potential 

This phenomenological potential was proposed by Song and Lin [16] 

VsL{r) = + Br'/^ + Co, (36) 
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with A = 0.923 GeV^/^ B = 0.511 GeV^/^ Co = -0.760 GeV, rUc = 1.820 GeV, and 
rrib — 5.190 GeV. The characteristic features of this potential may be traced in Ref. [16]. 



3. Turin potential 

Lichtenberg et al. [17] suggested such a potential which is an intermediate between the 
Cornell and Song- Lin potentials. It has the simple form 

Vrir) = -Ar~^/^ + Br^'^ + Co, (37) 

with A = 0.620 GeV^/^ B = 0.304 GeV^/^ Co = -0.783 GeV, rUc = 1.790 GeV and = 
5.166 GeV. The characteristic feature of this potential can be traced in Ref. [17]. 

4- Power-law potential 

a. Power-law potential of Martin [19,22].- The phenomenological power-law potential 

Vuir) = -8.056 GeV + (6.898 GeV)(r x 1 GeV)°-\ (38) 

is labeled as Martin's potential [19] with rUc — 1.80 GeV and rub — 5.162 GeV. 

b. Power-law potential of Rosner et al. [20.22]: 

Vti{r) = -0.726 GeV + (0.801 GeV) [(r x 1 GeV)^ - 1] /7, (39) 
with 7 = -0.12, mc = 1.56 GeV and = 4.96 GeV. 

5. Logarithmic potential of Quigg and Rosner [21,22] 

A Martin's power-law potential reduces [19] into 

VL{r) = -0.6161 GeV + (0.733 GeV) ln(r x 1 GeV), (40) 

with rUc — 1.50 GeV and = 4.890 GeV. The potential forms in (38), and (40) were used 
by Eichten and Quigg [1,21]. Further, all of these potential forms were also used for tjj and 
T data probing 0.1 fm < r < 1 fm region [15]. 
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B. QCD-motivated potentials 



1. Igi-Ono potential 

The Cornell, logarithmic and power law potentials for some systems should account 
properly for the running ag at the different distances. The and the string constant k 
that describes the long-distance behavior in (35) can vary independently. The Aqcd used 
to evaluate ag is related to a specific renormalization scheme so that comparison with other 
calculations is possible. This leads to expanding as to at least 2-loop order and use of 
nontrivial methods in the interpolation between long-and short-distance behaviors. The 
interquark potential at short distances has been calculated to 2-loop calculations in the 



modified minimal-subtraction (MS) scheme [25]. Together with the 2-loop expression for 
CKs, one has [10,25] 



V{r) 



3r 



27r 



{holniir + A) 



(41) 



and 



An 



bi In ln{jj,^/A- 



MS' 



(42) 



with 



[^'o, h, A] = 



2 38 31 5 ' 

11 - 102 - yn/, 6o7b + y ~ g^/ 



(43) 



Here n/ is the number of flavors with mass below and ■Je = 0.5772 is the Euler's number. 
The renormalization scale /i is usually chosen to be 1/r to obtain a simple form for V{r). 
The singularity in as{l/r) at r = l/Aj^ = 5 GeV"^ 1 f m for Aj^ = 200 MeV. This 
singularity can be removed by the substitution 



In 



r2A?_ 

MS 



f{r) = In 



r'^A- 

MS 



(44) 



The constant bo is an adjustable parameter of the potential and will not affect the pertur- 
bative part of the potential. Hence, setting n/ = 4 and n/ = 5 in Eq. (43), the one-gloun 
exchange part of the interquark potential simply take the following two forms 
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(ny=4) 



OGE 



(r) = 



and 



(«/=5), 



OGE 



r) 



167r 1 



IGtt 1 
' 23 r/(r) 



53 
75 



625 /(r) 



(45) 



^ _ 348 Zn/(r) ^27^ + 1 



(46) 



529 /(r) /(r) 

respectively, where /(r) is the function given in Eq. (44). Further, the long distance 
interquark potential grows linearly leading to confinement as 



VL{r) = oor. 



Igi and Ono [25] proposed a potential whose general form 



(47) 



^(n,=4,5)(^) ^ ^OG^"''^ + «or + dore-^'' + Co, 



(48) 



so as to interpolate smoothly between the two parts. They added phenomenologically a 
term dore''^°^ to the potential so that to adjust the intermediate range behavior by which 
the range of Aj^g is extended keeping linearly rising confining potential. Hence, the Ajj^ = 
(100, 500) MeV keeping a pretty good fit to the cc and bb data. 



2. Improved Chen-Kuang potential 

Chen and Kuang [26] proposed two improved potential models so that the parameters in 
Ref . [26] all vary explicitly with Aj^ so that these parameters can only be given numerically 
for several values of A^^^. Such potentials have the natural QCD interpretation and explicit 
^Ws dependence both for giving clear link between QCD and experiment and for convenience 
in practical calculation for a given value of Aj^. This potential takes the general form 



V^("/-)(.)= l^TT 1 



25 r/(r) 



53 -\ 



^ A62 Infjr) ^ 2je + % 



+ Kr, (49) 



625 /(r) /(r) 

where the string tension is related to Regge slope by k — The function /(r) in Eq. 
(49), takes the general form 
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/(r) = In 



-\ 2 



A 



+ 4.62 - A{r) 



MS' 



and 



A{r) 



1 - - 



lA 



MS 



4A:?— 
MS 



1 — exp 
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3^ - 1 



A 



MS' 



A 



MS' 



with the following set of parameters fit 



[k, a, A^^] = [0.1491 GeV^ 1.067 GeV'^, 180 MeV] . 



(50) 



(51) 



(52) 



together with rric = 1.478 GeV, = 4.876 GeV, Aj^ = (100, 600) MeV and = 0.235. 
The details of this potential can be traced in Ref . [26] . 



IV. SUMMARY AND CONCLUSIONS 

We have solved numerically SS wave equation, Eq. (3), for various potentials to determine 
the position of the spin-averaged ground-state masses of Mis(J/V' and r]c) — 3068 MeV, 
Mis(T and r^b) = 9447 MeV, and also Mis (5c) with the help of Eq. (29). For simphcity 
and for the sake of comparison, we have neglected the variation of as with momentum in 
(32) to have a common spectra for all states and scale the splitting of c6 and bb from the 
charmonium value. In our calculations for the hyperfine nsplitting values of the cb spectrum, 
we have followed the same approach of Ref. [9] in fitting the QCD coupling constant a^. 

Table I reports the SS SAD mass spectrum and hyperfine sphtting values of the vector 
and pseudoscalar masses of the cb for a family of seven central potentials. Masses and 
splittings lie within the ranges quoted by Kwong and Rosner [7], Eichten and Quigg [1] and 
also Fulcher et al. [4]. The value of the coupling constant ccs, is obtained from the proper 
splitting values of the cc and bb quarkonium systems [1,4,6-9]. 

It has been stated the possibility of producing c6 mesons in e"'"e~ and hadron-hadron 
colliders [31-34]. We have used the couphng constant in our analysis in the range 0.185 < 
as < 0.293 for all central potentials and 0.190 < < 0.448 for the QCD-motivated poten- 
tials as shown in Tables I and IV. Here, we have followed most authors in fixing the coupling 
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constant as{mc) in reproducing all the spectra, [c.f. e.g. Refs. [1,4,8,9]), although slight 
changes are made by other authors [15,35].^ 

Our results for SAD cb masses together with their lowest S-state splittings for both 
Schrodinger and Salpeter models are presented in Table II and compared with the other 
estimates of Refs. [4,24,32,36]. The values of the fits are hsted in Table II, they show 
that the Salpeter results yield = 447 in the present model and give better agreement with 
experiment, PDG [30], than the Schrodinger results of — 836. Thus, the of the overall 
fit increases by about a factor 2. Hence, we conclude that extending the contest between 
the Salpeter equation and Schrodinger's [9] to the charmonium and bottomonium systems 
provides additional evidence of an experimental signature for the semi-relativistic kinetic 
energy correction operator (second term) of Eq. (3). We have found that the measured 
values for the SAD of charmonium and the upsilon system are in better agreement with 
the Salpeter equation results than the Schrodinger equation results, (cf. e.g. Ref. [9]) and 
Table II of the present work, which is the same conclusion reached by Jacobs et al. [37] 
and Fulcher et al. [4]. Therefore, minimizing the quantity x^ is a stringent argument that 
indicates the accuracy of the semirelativistic model over the previous model [9] . Overmore, 
in Table II, we have considered two different fitted sets of parameters for reproducing the 
d) masses. We noted that the values of quark masses increase if we are allowed to add an 
additional constant Cq as given in the potential generality (33) and decrease if we drop it 
out. Hence, the overall fit increases by nearly a factor 4 in this trend. 

Bambilla and Vairo [3] have calculated the maximum final result of {Mb*)^^^^ = Q326tf 
MeV, the upper limit corresponds to the choice of parameters ^j^^ — 350 MeV and 
H = 1.2 GeV, while the lower limit to = 250 MeV and fj, = 2.0 GeV as the best 

approximation to their perturbative calculation. Our predictions for the cb fine and hy- 



^Kiselev et al. [15] have taken into account that AMt(IS) = ^^AM^(IS) with as{T)/as{tp) ~ 
3/4. On the other hand, Motyka and Zalewiski [35] also found ~ 11/18. 
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perfine splittings together with the ones estimated by other authors are hsted in Table III. 
The SLNET fine and hyperfine splitting estimations for a group potentials are all fall in 
the range demonstrated by other authors. Larger discrepancies among the various methods 
occur for the ground and excited states [6] . 

The fitted set of parameters and SAD mass spectrum and their splittings of the Igi-Ono 
potential [25] labeled by type 1, 2 and 3, are listed in Table IV, V and VI. It is clear that the 
overall study seems likely to be pretty good and the reproduced fine and hyperfine splittings 
of the states are also reasonable. Further, it is worthwhile to note that the quark masses rric 
and rub increase as the values of Aj^g increase. This is explained in the following way. The 
absolute value of the short-range asymptotic behavior of the potential in Eqs. (45) and (46) 
decreases with increasing Aj^. In order to reproduce the SAD masses we need larger quark 
masses for larger values of A^^^, (cf. e.g. Table III of Ref. [9]). Clearly, this is well noted in 
the rric values. 

We have also calculated the predicted cb spectrum obtained from the Chen-Kuang poten- 
tial in Table VII. We see that for states below the threshold, the deviations of the predicted 
spectra from the experimental SAD are within several MeV. We find that rric a-nd rrib are 
insensitive to the variation of Aj^ for this potential. The A;^ dependence of the Chen- 
Kuang potential is given in Eq. (49). The potential found to be more sensetive especially 
for the lowest states and found not sensitive for higher states. The effect of Aj^ is clearly 
on the Coulombic part of Eq. (49) . 

In Table VIII, in order to get some idea of an error estimate (precision) for our model, we 
have calculated the low-lying two S-state masses of the singlet and triplet Be system which 
are completely dependent on the running coupling constant ag values, (cf. e.g. Table I), 
determined in charmonium and the upsilon systems. Using the largest and smallest values of 
as in Table I to determine the errors, we have calculated Mb, = 6234^^^ MeV, Mb* = 6310^^ 
MeV, and 69 MeV < Ais < 80 MeV. Our hyperfine splitting value is in exellent agreement 
with the other estimates made by other authors. Clearly, the precision of the experiments [2] 
requires a very substantial improvement to be sensitive to the bound-state mass differences 
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between the various calculations listed in Table VIII. 

We study some properties of the system like the pseudoscalar decay constant given 
by the Van Royen-Weisskopf formula modified for color [4,38], that is. 



/Bc = i/^|i?is(0)|, (53) 

with the nonrelativistic radial wavefunction at the origin [21] 

|i?is(0)| = V4^|V^is(0)|. (54) 

Moreover, in Eq. (53), the Salpeter bound-state mass of the low-lying {n — 1,1 = 0) 
pseudoscalar Sg-state can be found via 

Mb, (0-) = me + mfc + 2Ei^o - 3 A£;hf/4, (55) 

and also the vector S*-state via 

Mb. (1-) = m, + + 2Ei,o + ^E^f/^- (56) 

whereas the square-mass difference can be simply found via 

AM^ = M|j (1-) - M|^(0-) = 2A£;hf [m^ + + 2E^,q - A£;hf/4] , . (57) 

with 

A£;hp = ^ |i?ls(0)|^ (58) 

Using our estimates of Mb, and Mb* in Table II, we find decay constant for SS and for S 
models. They are listed in Table IX together with the results of other authors [1,4,39-41]. 
All of these results are also in reasonable agreement with the lattice results [5] . 

The empirical result obtained by Collins et al. [42] for potential model wave functions 
at the origin, that is, 

\RbM? ^\Rj/mt'\Rrm'-\ (59) 
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provides another touchstone for our numerical work. Prom charmonium and upsilon calcu- 
lations, we get wave function at the origin for SS and for S as reported in Table IX which 
are in excellent agreement with the result of Fulcher [4], in which |i?s^(0)|^ ~ 1.81 GeV^. 
Moreover, the empirical relationship for the ground state hyperfine splittings, that is, 

- Mb^ ~ 0.7 {Mj/^ - M^f-'' (Mr - M^J°■=^^ (60) 

yields a sphtting of 63 MeV, about 1% and 5% lower than our estimations for S and SS 
equations, respectively, as shown in Table II. Both of these results are in reasonable agree- 
ment with the results of Collins et al. [42], Fulcher [4] and FCY [4] and other authors [39-41] 
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APPENDIX A: SOME USEFUL EXPRESSIONS 



Here we list the analytic expressions of P^^\ /3^^\ej, and Sj for the spinless Salpeter 



equation 




Y^^—^ + [(1 + 2n,) £2 + 3(1 + 2nr + 2nl)s^] 

[ej + 6(1 + 2nr) SiS^ + (11 + 30n^ + SOn^) e^] ^ 



— u 



(Al) 



(^m ^ [(1 + 2nr% + 3(1 + 2nr + 271^)^4 + 5(3 + 8n, + + 4n^)4 
- a;"^(l + 2nr)el + 12(1 + 2n^ + 2n^)£2£4 + 2 £1^1 
+ 2(21 + 59n^ + hlnl + 34n^)£^ + 6(1 + 2nr)eih 
+ 30(1 + 2nr + 2n^)£i4 + 2(11 + 30n^ + 30n^)£34 

+ 10(13 + AQur + 42n^ + 28n^)£3^5 + 6(1 + 2nr)eA] 
+ uj-^ [4 £^ £2 + 36(1 + 2n^)£i£2£3 + 8(11 + 30n^ + 30n^)£2£3 
+ 24(1 + 2nr)el£4^ + 8(31 + 78n^ + 78n^)£i£3£4 

+ 12 (57 + 189n^ + 225n^ + Ib^nDelsi] 

e\e^ + 108(1 + 2nr)elel + 48(11 + 30n^ + 30n^)£i£^ 
+ 30(31 + 109n, + 141n2 + 94n^)£|] , (A2) 



where 



(4//a;)V2' 



i = 1,2,3,4. 



(A3) 



and 



(4/xa;)j/2' 



j = l,2,3,4,5,6. 



(A4) 



(2-a) 
4// 



£2 = - — 




(A5) 




2y(ro)r"(ro) 
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GV'{ro)V"{ro) , 2V"'{ro)E^ 



m 



+ 



m 



(A6) 



^4 = T7^ \- 



16/1 24Q 



V""{ro) - 



2V{ro)V""{ro) 



m 



SV'{ro)V"'{ro) 6V" {ro)V" (vo) 2V""{ro)E, 



m 



m 



+ 



m 



(A7) 



5i 



(l-a)(3-a) , 2rlE2V'{ro) 



8/1 



+ 



m'Q 



(A8) 



So = 



3(l-a)(3-a) , ^E^V'iro) 



16// 



+ 



m Q 



(A9) 



(2-a) 5(2-a) 
(^3 = — ; ()4 = - 



2/x 



(AlO) 



S5 = + 



8// 120g 



V""'{ro) - 



2V{ro)V""'iro) 



m 



10V'{ro)V""{ro) 20V" {ro)V"' (ro) , 2V""'{ro)Eo 



m 



m 



+ 



m 



(All) 



7 



16// 720Q 



l^"""(ro) - 



2V{ro)V"""{ro) 12V'{ro)V""'{ro) 



m 



m 



2W"{r^)V"'{ro) ?,W"{To)V""{ro) 2V"""{to)Eo 



+ 



m 



m 



m 



(A12) 
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TABLE I. cb quarkonium masses and hyperfine splittings A„5*(in MeV) 
of the lowest two S-states calculated for a group of static potentials (GSP). 



States 


[6] 


Cornell 1 Cornell 2 Song- 


Lin Turin Martin Rosner Logarithmic 






0.293 


0.258 


0.242 


0.263 


0.235 


0.185 


0.200 


M{1S) 




6289 


6296 


6288 


6291 


6286 


6294 


6294 


M{l^Si) 


6327 


6312 


6313 


6308 


6311 


6304 


6312 


6311 






6220 


6244 


6230 


6231 


6232 


6238 


6241 




77 


91.7 


69.1 


77.9 


79.6 


72.2 


73.7 


69.4 


M{2S) 




6877 


6876 


6853 


6859 


6865 


6882 


6870 


M{2^Si) 


6906 


6886 


6885 


6861 


6868 


6874 


6889 


6878 






6850 


6848 


6827 


6832 


6836 


6859 


6846 




42 


36.1 


37.0 


34.1 


35.8 


37.6 


29.7 


32.4 


M{3S) 




7247 


7267 


7181 


7214 


7202 


7185 


7190 


M{4S) 




7545 


7586 


7422 


7491 


7444 


7387 


7412 


M(1P) 


6754 


6759 


6741 


6731 


6735 


6721 


6770 


6744 


M(2P) 


7154 


7133 


7140 


7089 


7109 


7102 


7113 


7104 


M(3P) 




7395 


7470 


7347 


7398 


7367 


7334 


7347 


M{1D) 


7028 


7018 


7014 


6998 


7002 


7003 


7043 


7021 


M{2D) 


7367 


7332 


7353 


7272 


7304 


7290 


7282 


7284 


M{3D) 




7562 


7649 


7488 


7556 


7509 


7454 


7479 



*Ans^M{n^Si)-M(n'So). 
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TABLE II. cb masses and A^^ (in MeV) of the using Coulomb- 
plus-linear potential (CPLP), in the Schrodinger and Salpeter models. 



States 


CPLP" 


CPLP*" CPLP*" [4]" [4]^ [32] [36] [24] 


asirric) 


0.276 


0.257 


0.322 




M(1S) 


6339 


6303 


6297 


6338 6314 6301 6317 6315 


M(l^Si) 


6355 


6319 


6317 


6332 




6291 


6253 


6235 


6270 




63.7 


66.2 


82.2 


62 


M(2S) 


6930 


6885 


6852 


6918 6872 6893 6870 7009 


M(2^Si) 


6940 


6894 


6862 


6881 


M(2^So) 


6902 


6856 


6823 


6835 


A25 


38 


38.5 


39.6 


46 


MiZS) 


7352 


7285 


7220 


7336 7270 7237 7225 


MiAS) 


7707 


7615 


7520 




M(1P) 


6756 


6741 


6725 


6755 6745 6728 6735 


M(2P) 


7195 


7152 


7100 


7193 7157 7122 


M(3P) 


7563 


7493 


7410 




MilD) 


7036 


7019 


6980 


7145 


M{2D) 


7418 


7371 


7299 




M(3D) 


7755 


7677 


7577 




Parameters 












U.4 ( z 


0.437 


0.457 






0.191 


0.203 


0.182 




CoiGeV) 


0.0 


0.0 


-0.790 




mdGeV) 


1.3205 


1.321 


1.79 




mb{GeV) 


4.7485 


4.731 


5.17 






836^ 


447c 


201V 





"In the Schrodinger model, ''In the Salpeter model, '^Fitting to cc and bb spectra. 
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TABLE III. Fine and hyperfine splittings of cb states (in MeV) using 
GSP, lOP and CKP in our work compared with the other authors. 
States Latt.[6] Qua. [6] Lin. [6] [4] lOP* GSP* CKP* 

Fine SpUtting 
M{2S) - M{1S) 
M{3S) - M{1S) 
M{4S) - M{1S) 
M{2P) - M{1P) 
M{3P) - M{1P) 
M(2D) - M{1D) 
M{3D) - M{1D) 
Hyperfine Sphtting 
Ais 

A25 

* The averaged sphttings in our work for different potentials. 



TABLE IV. Parameters used for Igi-Ono potential (lOP). 



ao (GeV') 


^0 (GeV) do (GeV') 


as{mc) 


rUc (GeV) 


rub (GeV) 


Type 1 (n/ 


= 5) 










0.1587 


0.3436 


0.2550 


0.190 


1.38 


4.782 


Type 2 (n/ 


= 4) 










0.1587 


0.3436 


0.2550 


0.403 


1.353 


4.746 


Type 3 (n/ 


= 4) 










0.1585 








0.448 


1.494 


4.885 



672 ± 120 558 533 579 577 ± 25 5781^^ 4941^9 

931 899 916+11 92lt|° 866l^| 

11761^? mstll^ lUdtll 

382 376 375 ±8 370^^^ 389 ± 1 

6571?^ GSltrl 6951? 

324 321 2961?^ 2881^^ 318^? 

547111 5Ut\ll 586 ±1 



41 ±20 77 62 55 63 ± 1 76^^^ 651^3 
30 ±8 42 33 32 32 35l^ 351?;^ 
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TABLE V. cb quarkonium mass spectrum and A„5 (in MeV) 
predicted by using lOP of types 1, 2, and 3 in our work. 



States Aj^{MeV) 


= 200 


200 


200 


250 


300 


400 


b 


20'' 


20^ 


23.3^ 


23.3" 23.3" 


23.3" 


M{1S) 


6295 


6294 


6364 


6324 


6292 


6245 


M{l^Si) 


6311 


6310 


6377 


6340 


6310 


6267 


M{l'So) 


6247 


6247 


6324 


6277 


6239 


6179 




64.6 


63.5 


53.3 


62.1 


70.9 


87.9 


M{2S) 


6897 


6872 


6899 


6876 


6857 


6829 


M(2^Si) 


6905 


6880 


6907 


6884 


6865 


6837 




6873 


6848 


6876 


6852 


6832 


6803 


A25 


32.2 


32.2 


31.2 


32.3 


33.1 


34 


M{3S) 


7231 


7197 


7249 


7233 


7222 


7210 


M{4S) 


7474 


7439 


7539 


7529 


7523 


7516 


M(1P) 


6757 


6733 


6755 


6741 


6731 


6719 


M(2P) 


7132 


7099 


7134 


7124 


7117 


7108 


M(3P) 


7404 


7368 


7437 


7430 


7425 


7420 


M{1D) 


7041 


7008 


7019 


7013 


7009 


7005 


M{2D) 


7328 


7292 


7335 


7330 


7327 


7323 


M{3D) 


7570 


7532 


7606 


7602 


7599 


7597 



" lype 1 with Co = -28 MeV to fit (c.o.g.) value. 

Type 2 with Co = -28 MeV to fit (c.o.g.) value 
" Type 3 with Co - -29 MeV to fit (c.o.g.) value. 
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TABLE VI. cb quarkonium mass spectrum 
and AM^ (in MeV) predicted by using lOP. 



States 


Ajg5(Mey) 200 


250 


300 


M(1S) 


6382 


6342 


6311 


M(2S) 


6918 


6894 


6875 


M(3S) 


7267 


7251 


7241 


M(4S) 


7558 


7548 


7542 


M(1P) 


6773 


6759 


6750 


M(2P) 


7153 


7142 


7135 


M(3P) 


7456 


7448 


7443 


M(1D) 


7037 


7031 


7027 


M(2D) 


7353 


7348 


7345 


M(3D) 


7625 


7620 


7618 


AM2S 


585^^ 536 


552 


564 


AMgs 


349 


357 


366 


AM4S 


291 


297 


301 


AM2P 


380 


383 


385 


AMgp 


303 


306 


308 


AM2D 


316 


317 


318 


AM3D 


272 


272 


273 



"Fine splitting is AM(„+i)s = M{{n + 1)S) - M(nS). 

''Type 3 fitted to the experimental cc and bb spectra with = 4.874 GeV. 
^Here we cite GKLT in Ref. [32]. 
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TABLE VII. cb quarkonium mass spectrum and A„5 



(in MeV) predicted by using CKP in 


our work. 








States 


Aj^^(MeV) 100 - 


270 275 - 


350 400 


450 


500 


600 


M(1S) 


6250 


6308 


6348 


6278 


6318 


6277 




6267 


6325 


6365 


6293 


6331 


6293 


M(liSo) 


6197 


6255 


6296 


6232 


6279 


6231 


Ais 


69.9 


69.9 


68.9 


68.9 


51.8 


61.6 


M(2S) 


6808 


6808 


6808 


6808 


6841 


6835 


M(23Si) 


6817 


6817 


6817 


6817 


6850 


6844 




6782 


6782 


6782 


6782 


6816 


6809 


Ass 


35.2 


35.2 


35.2 


35.2 


33.8 


35.1 


M(3S) 


7171 


7171 


7171 


7171 


7229 


7173 


M(4S) 


7466 


7466 


7466 


7466 


7490 


7466 


M(1P) 


6674 


6674 


6674 


6674 


6690 


6678 


M(2P) 


7062 


7062 


7062 


7062 


7080 


7065 


M(3P) 


7368 


7368 


7368 


7368 


7388 


7368 


M(1D) 


6953 


6953 


6953 


6953 


6957 


6955 


M(2D) 


7270 


7270 


7270 


7270 


7277 


7270 


M(3D) 


7540 


7540 


7540 


7540 


7542 


7540 



" Type 1 with Cq = —28 MeV to agree with (c.o.g.) value. 
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TABLE VIII. The predicted cb masses of the lowest S-wave and its 
sphtting (all in MeV) compared with the other authors. 



work 


M(1S) 




Mbj(135i) 


Ais 


Eichten et al. [ll 

L J 




6258 ± 20 






Colangelo et al. [3] 




6280 


6350 


70 


Baker et al. [33] 




6287 


6372 


85 


Roncafflia et al. [33] 






6320 ± 10 




Godfrey et al. fll 

o L J 




6270 


6340 




Bagan et al. [1,33J 




6255 ± 20 


6330 ± 20 


75 


Bambilla et al. [3] 






6326lf 




Baldicchi et al. [6] 




6194 ~ 6292 


6284 - 6357 65 < Ais < 90 


SLNET (GSP)t 


6291 ± 5 


62341^^ 


63101^ 


69 < Ais < 80 


SLNET (lOP)* 


6324 


6277 


6340 


62 


SLNET 


6291T 


6234 


6310 


66 



* Experimental mass of such a singlet state is presented in [2,6]. 

^ Averaging over the seven values in Table I. 

^ An estimation using type 3 with Ajg^ = 250 MeV. 

"I" Best estimation to the (c.o.g.) lowest S-state. 

TABLE IX. The leptonic Be meson constant and radial wave function 
at the origin calculated in our model and by the other authors. 

Level (l^^o) Martin Coulomb EQ[1] F[4] [39] [40] [41] S* SS* 
Jb, (MeV) 510 ± 80 456 ± 70 495 517 410 600 500 495 524 

|i?Be(0)|^ (GeV^) 1.716 - 1.638 1.81 - - - 1.85 1.90 

* Present work. 
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